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Abstract — In this paper, we study a discrete multispecies 
Lotka-Volterra competition system. Assume that the coefficients 
in the system are almost periodic sequences, we obtain the 
sufficient conditions for the existence of a unique almost 
periodic solution which is uniformly asymptotically stable by 
constructing a suitable Liapunov function. One example 
together with numerical simulation indicates the feasibility of 
the main results. 

Index Terms —Almost periodic solution, Discrete, 

Lotka-Volterra competition system, Permanence, Uniformly 
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I. Introduction 


In paper [1], Chen and Wu had investigated the dynamic 
behavior of the following discrete n-species Gilpin-Ayala 
competition model 

Xi(k+1) =Zi(fc)exp (1.1) 

t j= i J 


where i = 1, 2, • • • , n; x t (k) is the density of competition 
species i at k -th generation, a^k) measures the intensity of 
intraspecific competition or interspecific action of 
competition species, respectively, b^k) represents the 
intrinsic growth rate of the competition species x t . 6 are 
positive constants, b^k), a^k), /, j = 1, 2, • • • , n are all 
positive sequences bounded above and below by positive 
constants. Obviously, when 0^= 1, system (1.1) reduces to 
the traditional discrete multispecies Lotka-Volterra 
competition model 


Xi(k+1) = Xi(k)exp lbi{k) -^dijityxjik) (1.2) 

t 1 ' 

For general non-autonomous case, sufficient conditions 
which ensure the permanence and the global stability of 
system (1.1) and (1.2) are obtained; For periodic case, 
sufficient conditions which ensure the existence of an unique 
globally stable positive periodic solution of system (1.1) and 
(1.2) are obtained. 

Notice that the investigation of almost periodic solutions 
for difference equations is one of most important topics in the 
qualitative theory of difference equations due to its 
applications in biology, ecology, neural network, and so 
forth(see [2-13] and the references cited therein). Wang and 
Liu [3] studied a discrete Lotka-Volterra competitive system 




Xi (n + 1) = Xi ( n) exp 


**i(«) -ai( n )*i( n ) 


c 2 {n)x 2 (n) 
1 + x 2 (n) 


x 2 (n + 1) = x 2 (n) exp 


r 2 (n) - a 2 (n)x 2 (n) - 


ci(n)xi(n) 

1 + Xi(ll) 


n = 0.1,2, 


With the help of the methods of the Lyapunov function, 
some analysis techniques, and preliminary lemmas, they 
establish a criterion for the existence, uniqueness, and 
uniformly asymptotic stability of positive almost periodic 
solution of the system. However, few work has been done 
previously on an almost periodic version which is 
corresponding to system (1.2). Then, we will further 
investigate the global stability of almost periodic solution of 
system (1.2). 

Denote as Z and Z + the set of integers and the set of 
nonnegative integers, respectively. For any bounded sequence 
{g(ft)} defined on Z, define 

g u = sup g(n),g l = inf g{n). 

n£Z n ^ z 

Throughout this paper, we assume that: 

(HI) ajj(k ) and biik) are bounded positive almost 
periodic sequences such that 

0 < a ij < a ij{k) < a-, 0 < b[ < bi(k) < bf , i,j = 1,2,-** , n. 

From the point of view of biology, in the sequel, we assume 
that x(0) = (xi(0), x 2 (0), • • • , x n (0)) > 0. Then it is easy to see 
that, for given x(0) > 0, the system (1.1) has a positive 
sequence solution x(k) = (x\(k), x 2 (k), • • • , x n (k))(k ^ Z + ) 
passing through x(0). 

The remaining part of this paper is organized as follows: In 
Section 2, we will introduce some definitions and several 
useful lemmas. In Section 3, by applying the theory of 
difference inequality, we present the permanence results for 
system (1.2). In Section 4, we establish the sufficient 
conditions for the existence of a unique uniformly 
asymptotically stable almost periodic solution of system (1.2). 
The main results are illustrated by an example with a 
numerical simulation in the last section. 

II. Preliminaries 

First, we give the definitions of the terminologies involved. 

Definition 2.1([14]) A sequence x: Z^R is called an 
almost periodic sequence if the e-translation set of x 

E{e, z} = {t € Z :| x(n + r) — xf/z) |< e, Vn G Z} 

is a relatively dense set in Z for all s > 0; that is, for any given 
£ > 0, there exists an integer /(e) > 0 such that each interval of 
length /(e) contains an integer x^E{e, x} with 

x(n -hr) — x(n) \< e, Vn 6 Z. 

x is called an e-translation number of x(n). 

Lemma 2.1([15]) If {x(n)} is an almost periodic sequence, 
then {x(n )} is bounded. 

Lemma 2.2([16]) {x(n)} is an almost periodic sequence if 
and only if, for any sequence mf: Z, there exists a 
subsequence such that the sequence {x(n +m ik )} 
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converges uniformly for all n £= Z as k -*• 00 . Furthermore, the 
limit sequence is also an almost periodic sequence. 

Lemma 2.3([ 15]) Suppose that {pi(n)} and {p 2 (n)| are 
almost periodic real sequences. Then {pi(n)+p 2 (n)} and 
{pi(n)p 2 (n)} are almost periodic;l/pi(n) is also almost 
periodic provided that pi(n) ^0 for all n £ Z. 

Moreover, if e > 0 is an arbitrary real number, then there 
exists a relatively dense set that is s-almost periodic common 
to {pi(n)} and {p 2 (n)}. 

Lemma 2.4( [17]) Assume that sequence {x(n)} satisfies 
x(n) > 0 and 

x(n + 1) < x(n) cxp{a(n) — b(n)x(n)} 

for n^N, where a(n) and b(n) are non-negative sequences 
bounded above and below by positive constants. Then 

limsupx(fi) < yr exp{o Ti — 1}. 

n—H-oo O 

Lemma 2.5( [17]) Assume that sequence |x(n)} satisfies 

x(n -f- 1) > x(n) exp{a(n) — b(n)x(n)}, n > Nq, 


limsupx(n) < x* 

n— M-00 


and x(N 0 )>0, where a( n) and b(n) are non-negative sequences 
bounded above and below by positive constants and N 0 ^ N. 
Then 


lirninf x(n) > min 

n—> H-oc 




Consider the following almost periodic difference system: 

( 2 . 1 ) 

I x I < B}, and/(n, 


x(n + 1) = /(n, x(n)), n e Z + , 


where/: 77 xS B — R K , S B = {x <e R K : 
x) is almost periodic in n uniformly for x ^ S B and is 
continuous in x. The product system of (2.1) is the following 
system: 


x(n + 1) = f{n , x(n)), y{n + 1) = /(n, y(n)), 


( 2 . 2 ) 


and Zhang [18] obtained the following Theorem. 

Theorem 2.6( [18]) Suppose that there exists a Lyapunov 


x 


< B. 


function V (n, x, y) defined for n £ Z , 
satisfying the following conditions: 

(i) o(|| x-y ||) < V(n,x,y) < 6(|| x-y 
where a, b G K with K = {a G C(/? + , i? + ) : a(0) = 0 
and a is increasing}; 

(ii) || V(n,x u y x ) - V(n,x 2 ,y 2 ) ||< L 


y 


< B 


2*, - X n 


+ II Vi ~ y -2 II]? where L > 0 is a constant; 
(iii) AV( 2.2 ){ n , x ,y) ^ —q V (/?, x, y), where 0 < a < 1 
is a constant, and 


AV( 2 . 2 )(n,x,y) = V(n + 1 , f(n,x), f(n,y)) - V(n,x,y). 

Moreover, if there exists a solution cp(n) of (2.1) such that || 
cp(n) | ^ B*< B for n £ Z + , then there exists a unique 
uniformly asymptotically stable almost periodic solution p(n) 
of system (2.1) which is bounded by B*. In particular, if/(n, x) 
is periodic of period co, then there exists a unique uniformly 
asymptotically stable periodic solution of system (2.1) of 
period co. 


III. Permanence 

In this section, we establish a permanence result for system 
(1.2), which can be found by Lemma 2.4 and 2.5. 

Proposition 3.1 Assume that (HI) holds. Then any positive 


solution (xi(k), x 2 (k), • • • , x n (k)) of system (1.2) satisfies 
limsup Xi(k) < Mi, (3.1) 

Ic —^-j-OO 

where 




Proposition 3.2 Assume that (HI) and 

(H2) b\- ± a%jMj > 0 


1 


71. 


hold for all i = 1, 2, • • • , n, where M i? i = 1, 2, • • • , n are 
defined by (3.1). Then for every solution (x^k), x 2 (k), • • •, 
x n (k)) of system (1.1) satisfies 

Urn inf Xi(k) > rrii , 


where 


777 i = min 


n 


b\ - 


n 


E 

-exp 16' V 

*■ j=i 
n 


a u 

u lt 


b[- E afjMj 


3= 


a 


u 

ii 


i = 1,2, •• • ,n. 


Theorem 3.3 Assume that (HI) and (H2) hold, then system 
(1.1) is permanent. 


The next result tells us that there exist solutions of system 
(1.2) totally in the interval of Theorem 3.3. We denote by Q 
the set of all solutions (x^k), x 2 (k), • • • , x n (k)) of system (1.2) 
satisfying mi ^ Xi(k) ^ Mi(i =1, 2, • • • , n) for all k £ Z + . 
Proposition 3.4 Assume that (HI) and (H2) hold. Then 
O. 

Proof. By the almost periodicity of {ay(k)} and {bi(k)}, there 
exists an integer valued sequence {5 p } with 8 p ^+oo aS p-* 
+°° such that 

A* T Sp) —> ciij(k), bi(k + S p ) —> bi(k), as p —-Hoc. 

Let 8 be an arbitrary small positive number. It follows from 
Theorem 3.3 that there exists a positive integer N 0 such that 

nit — e < Xi(k) < Mi + £, k > No. 

Write x ip (k) = Xi(k + 5 P ) for k ^ N 0 - 5 P and p = 1, 2, • • • . 
For any positive integer q, it is easy to see that there exists a 
sequence {x ip (k) : p ^ q} such that the sequence x p (k) has a 
subsequence, denoted by |x ip (k)} again, converging on any 
finite interval of Z as p-* 00 . Thus we have a sequence {yi(k)} 
such that 

Xip(k) —> yi(k) for k G Z + as p —> -foe. 

This, combined with 
Xi(k + 1 + S p ) = Xi(k + S p ) exp j b t (k + S p ) 

n 

^ Sp)Xj (A “b Sp) 

3 =1 

gives us 

2/i(fc+ 1) = Vi{k) exp | bi(k) I- ^ Q,j(fc)yj(fc) Lz = 1,2, • • ,n. 

( j= i J 

We can easily see that (yi(k), y 2 (k), • • • , y n (k)) is a solution 
of system (1.2) and m { - 8 ^yi(k)^ Mi + 8 for k £ Z + . Since 8 
is an arbitrarily small positive number, it follows that m A ^ 
yi(k)^ Mi and hence we complete the proof. 
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IV. Almost Periodic Solution 
The main results of this paper concern the existence of a 


Moreover, for any 
(k,Q,W),(k,Q,W) eZ+ x S* x S*, 


unique uniformly asymptomatically stable almost periodic we have 


solution of system (1.2) by constructing a non-negative 
Lyapunov function. 

Theorem 4.1 Assume that (HI), (H2) and 

(H3) 0 </? = min {ft} < 1 

l<i<n 

hold, where 


\V(k,Q,W)-V(k,Q,W)\ 

n n 

J2(Pi( k ) ~ Qi( k )f ~ - m{k)) 


r=i 

n 


\PiK K ) ~ " n '" 2 

i=i 


< 1 (M k ) - Qi ( k )) 2 - ( Pi( k ) - Qi( k )) 


i= l 
n 


n 


Pi = 2a',m t - < t 2 M 2 - £ [M 2 a“? + (1 + 

n 

+ (1 + + MiMj Y. “«<*«] - 

i = 1, 2, • • • , n. Then there exists a unique uniformly 
asymptotically stable almost periodic solution (x^k), x 2 (k), 

• • • , x n (k)) of system (1.2) which is bounded by Q for all k e 
Z + . 

Proof. Let pi(k) = In Xi(k), i = 1, 2, • • • , n. From system (1.2), 
we have 

71 

Pi(k + 1 ) =Pi(k) +bi(k) - y^aij(k)e PiW , i = 1 , 2 ,- 


= I(PfW “ Qi( k )) + ( Pi ( k ) ~ Qi{ k ))\\(Pi( k ) - Qi( k )) - ( Pi( k ) ~ Qi ( k ))I 


i= l 
n 


< + h( k )\ + \pi( k )\ + l 9 iWI)(|p*W -PtWI + k«(fc) -g«(fc)l) 


i=i 


n 


< ^[51 \Pi( k ) ~ M k )\ +Y 1 M k ) 9 i( fc )l] 


i= 1 


i— 1 


= MIIQ - Q II + liw - m 

where 


Q = (pi(fc),p 2 (fc),-*- ,Pn(*0), IF = (^1 (A:),q 2 (^). 
and 

L = 4max{ max {A*}, max {IL}}. 

l<i<n l<i<n 


,q n (k)), 


, li. 


(4.1) 


3 =1 


From Proposition 3.4, we know that system (4.1) have 
bounded solution (pi(k), p 2 (k), • • • , p n (k)) satisfying 

In rrii < Pi(k) < In M*, i = 1,2,--- , n, k G Z~. 

Hence, lpi(k)l ^ A i? where Aj = max{I In mj, I In Mil}, i = 1, 
2, • • • , n. 


Thus, condition (ii) of Theorem 2.6 is satisfied. 

Finally, calculating the AV(k) of V(k) along the solutions 
of system (4.2), we have 

AV {4 . 2) (k) = V(k + l)-V(k) 

n n 

= ^2(Pii k + 1) - Qii k + l)) 2 - ^(Pi( k ) ~ Qi(k)) 2 


t=l 

n 


i=l 


= ^[{Pi{ k + 1 ) - Qi( k + l )) 2 - (P*(fc) “ 9 tW) 2 ] 


n 


For X ^ R n , we define the norm X = ^ x i 


z=l 

n 


n 


l 2 


i=l 


= X] 0>i(*)-«(*))-osi(*)(e p,( * ) -e , ‘ ( " !) )+ XI “<jW(e p,(fc) -e* (lb) ) 


1=1 


Consider the product system of system (4.1) 

n 

Pi(k + 1) = Pi{k) + bi(k) - ^2 a ij(k)e Pj ( k \ 

j =i 

n 

q t (k + 1 ) = qi(k) + b t (k ) - , i = 1 , 2 ,••• ,n. 

J=l 

We assume that Q = (pi(k), p 2 (k), • • • , p n (k)), W = (q,(k), 


- (Pi(fc) -Qi(k))' 


(4.2) 


= X 1 - 2au(k)(Pi(k) - < ?i (*))(e pi( *•' , “ + 4(k)(e p,(k> - e ,l(fc) )' 


i=l 


n 


— 2 a,i(A')(e p <( fc) E a ij(k)(e Pj ( fc ) — 

2 


” v anouiiiv Liiu-i, 'V — VF1V A V> j-/ 2 V-*V5 , ty — VH1V A V> / ” \ 

q 2 (k), • • • , q n (k)) are any two solutions of system (4.1) + y H a ti (/c)(e Pj(A ' ) - 

defined on Z + xS*; then, || Q || ^ B, W ^ B, where B ; 1 


n 


= X {+ Bi}, s* = {(Pi(k), p 2 (k), • • • , p n (k))l In mi^ Pi(n) 


+ 2 ( Pi (k) - qi (k)) X <Hj(*)(e»«-e®W)}. 


i=l 


^ lnMi,i= 1,2, • • • ,n,k G Z + }. 

Let us construct a Lyapunov function defined on Z + x S' 
x S* as follows: 


71 , 


n 


v(k.Q.w) = Y t (Pi(k)-q^)) 2 - 


By the mean value theorem, it derives that 

e P.(fc) _ e *W = ^(k)( Pi (k) - q t (k)), i = 1,2, -- 

where £,( !<) lies between e pi(k) and e qi<k) . Then, we have 

AV ( 4 . 2 )« = X ( - 2 au(k)ii(k)( Pi (k) - q,(k )) 2 + (k)(p,(k) - q t (k)f 


1=1 


1=1 


n 


n 


It is obvious that the norm || Q-W || = S lpi(k)—qi(k)l is 

i=l 


n 


- 2 au{k)£i(k)(pi(k) - q { (k )) a«j( fe Kj( fc )(Pj( fc ) “ 9 j( fc )) 

j=i.jy* 

+(ii 


equivalent to | Q-W *= £ (pi(k)-qi(k)) ] ; that is, there 

i=l 

are two constants Ci > 0, c 2 > 0, such that 

ci\\Q-W\\<\\Q-W\\*<c 2 \\Q-W\l 


+ 


2 {Pi{ k ) ~ Qi{ k )) ^ a ij( k )tj( k )(Pj( k ) ~ 9 j(^))| 
j=i.jy* 


= E 


n 


i=l 


2 a„(A’)C(^) + a%(k)g(k) + X a ?*(% 2 ( fc )) ~ 


then 


Let 


(ci||Q - IT||) 2 < T(A-, Q, W) < (c 2 \\Q - W\\f 


+ 2 X [[1 -<>«(*)?•(*)]“«(%(*) + ^ X “iiC*>€i(*)«ij(*)€j(*)) x 


(Pi( k ) - Qi(k)){Pj(k) - qj{k )) 


• 0 , G C(/? + , f? + ), - 0 (t) = t^(x) = c|t 2 ; 

then, condition (i) of Theorem 2.6 is satisfied. 
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n 


< Y \ ( - 2 a ii(k)^(k) + al{h)g{k) + Y a^k^ik) j ( Pi {k) - qi(k)) 2 


i= 1 


+2 


Y (l 1 - a ii( k )£i( k )] a ij( k )Zj( k ) + | Y a li( k )^i( k ) a u( k K( k )) X 
{Pii k ) - qi( k ))ipj{ k ) - Qj{ k )) }• 

Then, we have 

n 

AV ( 4 . 2 ) (fc)<^[V u (fc) + V i 2 (A-)], 

where 


i=l 


V i2 (k) = 2 


v n( k ) = ( - 2 aii(k)€i(k) + a? f (fc)^(fc) + ^ a 2 t {k)(?{k) j (. Pi (k ) - ?i(fc )) 2 

< ( - 2 < 4 m * + + M ? Y a j > 2 ) (&( k ) - &( k )) 2 ’ 

E ([l-a„(*)&((r)]ay(fcK Jf (A)+i £ <>«(*)&(*)««(*)«#)) x 

J=1.J7 4 * ' 1=1 .i#*J ' 

(Pi(/c) - qi(k))( Pj (k) - qj(k)) 

< E ((1 + a“ j M j )a“-M,+ iiW, M, E ««««) 1<P»W - «<*))* + <*»(*)- *(*))*) 

Hence, we have 

n / / n . vi v 

AV' ( 4.2)W < H (-2a' t m,+ a “ 2 A/ 2 + ^ [A/^ + (l+a5M 4 )a“ i A/ > + -M j A/ j ^ «,>,“]) x 

iP,(k) - q,(k)f + 53 ((l + a-A/.K^ + ^M/j 53 ofo-W*)-^*)) 2 } 

n ( / ”r l”\ 

= 53{ ( - 2a!,-", + a" 2 A/, 2 + 53 [A/ 2 a“ 2 + (1 + aJAf.KjA/j + — A/,A/j 53 a,“a ( “ ] ) (p,(Ar) - q,(k)f 

+ 53 ((1 + a^A/, )a^A/, + ^M,Mj 53 a,x)(P.(*)-V.(A-)) 2 } 

n / / n n \ 

= 53 ( - 2o " m ' + a “ 2 A/ 2 + 53 [Mfatf + (l + fl“ A/,)«“ A/j + (1 + + M.M, £ ) * 

i=l '' j -\JfH L l=U»Hj ' 

(Pi(*)-9i(*)) 2 | 

n / / n n v 

= - 53 ( 2a » m > - a “ 2jU - 2 - Z [ U " a “. 2 + 0 + aSA/<KA#i + (1 + a* j Af i )a'‘ I M i + A/.A/, 53 ajSafc] ) x 

(p.(fr) - q,(k)f | 

n 

<-53^(p ( (fc)- ft (jt)) a 

«—i 

ra 

< -^53(p<(*) ~ 9«(A - )) 2 
1=1 

= -/JV(*.Q,WO, 

where p = min {Pi}. That is, there exists a positive constant 0 

l<i<n 

< P < 1 such that 

AV(4.2 ){k,Q,W) < —f3V (k,Q, W). 

From 0<P <1, the condition (iii) of Theorem 2.6 is satisfied. 
So, according to Theorem 2.6, there exists a unique uniformly 
asymptotically stable almost periodic solution (pi(k), 
p 2 (k), • • • , p n (k)) of system (4.1) which is bounded by S* for 
all k ^ Z + . It means that there exists a unique uniformly 
asymptotically stable almost periodic solution (x^k), 
x 2 (k), • • • , x n (k)) of system (1.2) which is bounded by Q for 
all k ^ Z + . This completed the proof. 2 
Remark 4.2 If n = 2, the conditions of Theorem 4.1 can be 
simplified. Therefore, we have the following results. 
Corollary 4.3 Let n = 2, assume that (HI), (H2) and 

0 < f$ = min{/?i2,/?2i} < 1 

hold, where 

Pij = 2 Jum - afMf - Alfa]? - (1 + 2 

- (1 + 2ajj Mj )a“jA/j, 

i, j = 1, 2, j^i. Then there exists a unique uniformly 
asymptotically stable almost periodic solution (x^k), x 2 (k)) 
of system (1.2) which is bounded by Q for all k e Z + . 


(5.1) 


V . NUMERICAL SIMULATION 

In this section, we give the following examples to check the 
feasibility of our results. 

Example 5.1 Consider the discrete multispecies 
Lotka-Volterra competition system: 

X\(k -f 1) = xi(k)oxp 11.1 — 0.01 sin(\/3 k) — (1.15 — 0.01 sin(v'2A:)):ri(fc) 

-(0.055 + 0.(K)2cos(v / 5A - ))x 2 (^) - (0.03 + 0.002cos( v / 5Jt))3' 3 (Jt)|, 
or 2 {k + 1) = r 2 (A-)exp |l.1 — 0.025sin(y/3fc) — (0.02 — 0.003cos(v/5A-))T|(A*) 

-(1.08 + 0.025 sin{y/2k))x 2 (k) - (0.025 + 0.002ct*(v/5fc))z 3 (*) j, 
x 3 (k + 1) = x 3 (k )exp j 1.15 — 0.02sin(\/2A:) — (0.04 + 0.0025cos( \ r Ak))x l (k) 

- (0.029 + 0.0012 sin( y/2k))x 2 (A) — (1.13 -J- 0.02 sin( \/5n ) )x 3 (it) J. 

A computation shows that 

77ii « 0.7743, Mi « 0.9807, rn 2 ~ 0.9109, 

M 2 ~ 1.1572, 7713 ~ 0.8123, M 3 « 1.0579, 
and moreover, we have 

ft w 0.0652, ft ~ 0.0297, ft w 0.0451, 

that 0 < min{p 1? p 2 , p 3 } < 1. It is easy to see that the condition 
(H2) and (H3) are satisfied. Hence, there exists a unique 
uniformly asymptotically stable almost periodic solution of 
system (5.1). Our numerical simulations support our 
results(see Figs. 1,2 and 3). 



FIGURE 1: Dynamic behavior of the first component Xx(k) of 
the solution (x^k), x 2 (k), x 3 (k)) to system (5.1) with the initial 
conditions (0.87,1.02,1.03), (0.93,1.13,0.86) and 

(0.81,0.97,0.97) fork £ [1, 80], respectively. 

1.15 


X 



1.05- 


0.95 

0 10 20 30 40 50 60 70 80 

time k 

FIGURE2: Dynamic behavior of the second component x 2 (k) 
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of the solution (xi(k), x 2 (k), x 3 (k)) to system (5.1) with the 
initial conditions (0.87,1.02,1.03), (0.93,1.13,0.86) and 
(0.81,0.97,0.97) fork £ [1, 80], respectively. 



FIGURE3: Dynamic behavior of the third component x 3 (k) of 
the solution (x^k), x 2 (k), x 3 (k)) to system (5.1) with the initial 
conditions (0.87,1.02,1.03), (0.93,1.13,0.86) and 

(0.81,0.97,0.97) fork £ [1, 80], respectively. 
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